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ABSTRACT 

Context. We present the results of a near-infrared photometric and spectroscopic study of the star forming region G6 1.48+0.09. 
Aims. The purpose of this study is to characterize the stellar content of the cluster and to determine its distance, extinction, age and 
mass. 

Methods. The stellar population was studied by using color-magnitude diagrams to select twenty promising cluster members, for 
which follow up spectroscopy was done. The observed spectra allowed a spectral classification of the stars. 
Results. Two stars have emission lines, twelve are G-type stars, and six are late-0 or early-B stars. 

Conclusions. The cluster's extinction varies from Aks = 0.9 to A^s = 2.6, (or Ay ~ & to Ay ~ 23). G61.48+0.09 is a star forming 
region located at 2.5±0.4 Kpc. The cluster is younger than 10 Myr and has a minimum stellar mass of 1 .5±0.5x 10^^ Mq. However, the 
actual total mass of the cluster remains undetermined, as we cannot see its whole stellar content. 
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> 1. Introduction 

The exploration of the Milky Way is experiencing a strong 
' boost due to th e use of recent infrared surveys, like 2MAS S 
^ ; (ISkrutskieetalJ 12006,1 and GLIMPSE ( Benjamin etal J l2003h . 
■ frequently combined with data from other wavelength ranges, 
' like radio or y-rays. A strong effort is being dedicated to the dis- 
0^ covery of n ew stellar clusters , containing massive stars (Ar ches, 
O Quintuplet, iFiger et alJ (l99% . RSGCl, Figeretal.1 (1200^ 1. or 
to a deeper study of already known stellar clusters and cluster 
• candidates, resulting sometim es in their identifi cation as mas- 
sive stellar clusters (Cyg OB2. iKnodlsedeij ( 12000 '). Westerlundl, 
H LClark & Negueruela (2002)). The search has been st rongly en- 
- - ■ abled by new catalogs of st e llar cluster candidate s Uke lBica et al. 



may play a substantial role in the star formation history of the 
Milky Way and other galaxies. 

This is the first paper of our MASGOMAS (MAssive Stars 
in Galactic Obscured MAssive clusterS) project, a study of mas- 
sive star cluster candidates in the Milky Way that is being carried 
out using LIRIS at WHT. Our idea is to select known clusters 
or cluster candidates from the literature that could turn out to 
be massive clusters, observe them photometrically and perform 
a spectroscopic follow-up of the most promising ones. In this 
paper we present a photometric and spectroscopic study of the 
stellar po pulation of G61.48+0.09, a cluster that has got recent 
attention Pug a et al] (120041) and whose distance is matter of de- 
bate (see below). 

G61.48+0.09 is a complex star-forming region located in 
the direction of the Cygnus arm. Figure [T] shows a near- 
(left panel) and a mid- (right panel) infrared images of the 
G61.4 8-h0.09 region. Co ntinuum emission contour levels at 8.3 
GHz dGarav et alJll998i) are also shown. As it can be seen in 
the right panel of this figure, the G61.48+0.09 region consists 
of some bright stars surrounded by a HII region, Sh2-88, com- 
posed of two sub-regio ns, Bl (westwards) and B2 (eastwards; 
iDeharveng et alJ l2000l) . The central part is a triplet of stars. 
Although the region has been extensively studied by different 
authors, no spectroscopy of their stars exists. Puga et al. (200-^, 
in their polarimetric study of the region,present a detailed review 
of previous radio and near-IR photometric work. From their po- 
larimetric study and th e energetics balance they conclude (as did 
IDeharveng et alj QoOO*) that the main ionizing source of this re- 
gion is probably still hidden behind obscuring clouds of natal 
material. A strong source in the L' band, identified as LI, and a 
fainter L-source identified as L2, are proposed as possible main 



iBica et al 

(l2003h. iDutra et alj (l2003h . iMercer et al.l (|2005|) and iFroebrich 



( l2007h . This remarkable effort is leading us towards a re- 
evaluation of our ideas about the Milky Way star-forming ac- 
tivity (see f. e Figer 2008). 

iHanson & Popescul (l2008l) review the recent advances, and 
predict that more clusters are to b e discovered or ident ified. 
Recent exampl es, like Messineo 1 (iMessineo et al.l [20081) and 
DBS2003-179 dBorissova et al.ll2008l) . seem to confirm this pre- 
diction. A difference with former results, however, is that both 
clusters are moderate in their masses (about 2000 and 7000 M©, 
respectively), which indicates that we are discovering represen- 
tatives of a probably large number of intermediate mass stellar 
clusters. This is consistent with what we could expect, as the 
luminosity function of clusters in spiral galaxies follows a de- 
clining power law with exponent a= -2 (Larsen 2002), which 
for clusters of similar age can be translated into the same law for 
the mass function. These clusters, being massive and numerous. 
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ionizing sources. LI and L2 locations have also been marked in 
Figure [T| In the right panel, LI can be clearly seen as a bright 
object in the center of the Hll subregion Sh2-88B2. 

The distance to G61.48+0.09 is somewhat controversial. 
The first determ i nation was a kinemat ic distan c e deri ved by 
iReifenstein et all ( Il970l) . who using the ISchmidj ( 1 19651) model 
of the Galaxy arrived at the two possible distances of 2.0 and 7.5 
Kpc. Since then, many authors have made new determinations 
obtaining different values, but always in roug h agreement with 
one of both values obtained by Reifenstein. Churchwell et alJ 
(119901) give a distance of 5.4 kpc based on ammonia observa- 
tions combined with th e identifi cation of the region with CO ob- 
servations from Solomon et alJ (11987) . who use different tech- 
niques to discriminate clo se and far distances i n the list of gi- 
ant molecular clouds from ' Sanders et alj (Il986 i). A com parison 
of d ifferent determinations may be found in .Deharve ng et alJ 
( 1200 0). that adopt the short distanc e of 2.5 kpc, as most a uthors 
do. Nevertheless, the distance fr omlCh urchwe llet al.l ([1990) has 
very recently being adopted by IZhu e t al. (20Qa) in their study 
of Galactic clusters. Of course, many derived global properties 
of G6 1.48 -1-0. 09, like its total stellar mass or energetic content, 
may strongly depend on the adopted distance. 

We present new NIR photometric data and, for the first 
time, NIR spectroscopy of stars in the cluster region. Section |2] 
presents the observations and data reduction. In Section [3] we 
present the results of the photometric ( 13.1b and spectroscopic 
( I3.2l i analysis. Section |4]presents the discussion of these results 
and in Section|5]we present our conclusions. 



2. Observations and data reduction 

This study is based on broad-band imaging, long-slit and multi 
object spectroscopic observations done with LIRIS, a near-lR 
imager/spectrograph mounted at the Cassegrain focus of the 
4.2m William Herschel Telescope (La Palma). Table [T] gives a 
summary of the observation details. 

2.1. Photometry 

LlRlS is equipped with a Hawaii 1024 x 1024 HgCdTe aiTay 
detector, with a spatial scale of 0.25 "pixet^, providing a field 
of view of 4.27 ' on a side. Broad-band imaging observations 
were performed using the J, H and Ks filters. In order to opti- 
mize cosmic rays/bad pixels rejection and sky computation, a 
standard dithering mode (8-point pattern) was used. Two dif- 
ferent fields have been observed: the target field, centered in 
the G61.48+0.09 region, and a control field located ~8' west 
from the target field. Figure |2] shows a DSS-2-infrared survey 
image in which both target (a) and control (b) fields have been 
marked. It can be seen that there is a lack of stars around the 
G61.48+0.09 field if compared to the control filed. This sug- 
gests that there is a cloud associated to the G6 1.48 +0.09 region 
which is responsible for the increased extinction of the back- 
ground field stars. The data were reduced using the FATBOY 
(lEikenberrv et"ani2006t l^rner et al.ll2009 ') and LIRIS-QL im- 
age reduction packages following standard NIR reduction pro- 
cedures. The average seeing during the imaging observations of 
the G61.48+0.09 region was ~0.6". In the left panel ofFigurefT] 
the resultant G6 1.48 +0.09 image in the Ks -band is shown. 

DAOPHOT 11, ALLSTAR, and ALLFRAME ( Stetson el alj 
11994 ) were then used to obtain the instrumental photometry of 
the stars. The photometric lists of the observed targets were 
cleaned of non-stellar objects and poorly measured stars on 



the basis of the SHARP and PSF fitting cr parameters pro- 
vided by ALLFRAME. Only stars with |SHARP| < 0.25 and 
cr < 0. 1 have been considered. The photometric calibration was 
based on the 2MASS photometry. Figure |3] shows the 2MASS 
color-magnitude diagram (CMD) of the target field. For the cal- 
ibration, a set of isolated, non-saturated stars for which good 
photometry exists in both 2MASS and our images was selected. 
These stars have been marked in Figure|3]with open circles. They 
span a wide range of magnitudes (13 > /Ts > 9 mag) and col- 
ors {Q < { J - Ks) < 5). Comparing their instrumental photom- 
etry with the 2MASS catalogue, a calibration equation (which 
includes a color term) was derived for the /, H and Ks filters. 
Figure ID shows the photometric calibration results for our target 
field. Left panels show the photometric transformation from in- 
strumental (INSTR) to calibrated 2MASS magnitudes. Central 
panels show the difference between the calibrated photometry 
and the 2MASS photometry. It can be seen that this difference is 
smaller than 0.05 mag for most of the selected stars. Right pan- 
els show the PSF fitting cr parameter provided by ALLFRAME. 
This calibration procedure has been repeated for the control 
field. Finally, as the few stars brighter than Ks ~ 10 appear sat- 
urated in our data, their magnitudes have been obtained directly 
from the 2MASS catalogue. 

It is worth mentioning that the 2MASS photometry for most 
of the stars brighter than Ks ~ II and redder than (J - Ks) ~ 
is strongly affected by the molecular cloud. For those stars, the 
2MASS catalogue is adding the emission of dense cloud clumps 
to the star photometry, making the star artificially brighter and 
redder. In fact, some of the 2MASS detections in this region of 
the CMD are not real stars, but molecular cloud clumps instead. 
This effect is due to the limited 2MASS spatial resolution. For 
this reason, it is not possible to increase the number of stars with 
(J - Ks) >~ A--5 used for the photometric calibration. The cali- 
brated photometry is then less reliable for stars with ( J-Ks ) > 5. 
But, as it will be shown in Sections[3]and|4l these stars are back- 
ground late type stars that have not been used for the cluster 
description, with the exception of a YSO. As a conclusion, this 
lack of highly reddened stars in the photometric calibration is 
not affecting cluster results significantly. 

2.2. Spectroscopy 

We obtained medium-resolution HK- and /T-band spectra with 
WHT/LIRIS in multi-object (MOS) and long-slit (LS) spectro- 
scopic modes. The slit width was 0.8", allowing a spectral reso- 
lution of A/AA ~ 945 and 2500 in the HK- and /T-band spectra, 
respectively. We observed in a standard nodding pattern ABBA 
with individual exposure times of 20, 40 and 120 s (HK) and 
120 s (K). Table[T]lists the details of the spectroscopic observa- 
tions, and Fig. |5] shows a false color composition image of the 
G61.48+0.09 field marked with the position of the spectroscop- 
ically observed stars (see Section l372l l. 

The observations were reduced as follows: We first scaled B 
fields to their corresponding A fields using the quotient of me- 
dians as scaling factors. This step crudely corrects for average 
sky variations between frames with large integration times. We 
then subtracted B frames from A frames, and A frames from 
B frames before applying flatfield correction. The AB and BA 
frames were shifted in the spatial direction so that all positive 
continuum traces coincide before all frames were combined into 
a master frame. We then extracted the spectra from the master 
frame. The spectra were correlated with that of the standard star 
yielding wavelength solutions accounting for distortions of high 
orders. This solution was then applied to the spectra before divi- 



Marin-Franch et al.: The stellar population of the star forming region G6 1.48+0.09 



3 



sion by the spectrum of the standard star from which the stellar 
features had been removed. 



3. Results 

3.1. Color-Magnitude diagrams 

The CMDs of the target and control fields can provide valuable 
information about the G6 1.48 +0.09 region's extinction, distance 
and stellar content. Figure |6] shows the CMDs corresponding to 
the G61.48+0.09 region (a) and the control field (b). 

The position of a star on the CMD is determined from its ab- 
solute magnitude, intrinsic color, distance and extinction. When 
comparing target and control fields CMDs, the first major dif- 
ference that can be observed is the absence of the red-clump 
strip in the target field, while this feature can be clearly seen 
in the control one. The red-clump strip is the CMD trace that 
is populated by red-clump sta rs located at different distances. 
iLopez-Corredoira et al.l ( |2002|) developed a method to use red- 
clump stars as a distance indicator. In essence, this method uses 
a K2III star as representative red-clump star and, assuming an 
extinction law, computes the expected position of this star in the 
CMD when located at different distances. Comparing the for- 
mer with the observed red-clump strip, the distance of red-clump 
stars can be estimated. 

In this paper w e adopt the|Coi(|2000!) absolute V-magnitude 
and iDucati et all (|2001|) IR colors for a K2III star The black 
solid line in Figure |6] shows the expected position of a K2III star 
as a function of i ts distance. The Rieke et al. ( 1989) extinction 
law and R=3.09 ( IRieke & Lebofskvl il985) have been assumed. 
The fact that the red-clump strip is not seen in the G61 .48+0.09 
region shows that the cluster, in particular its associated cloud, 
is adding extinction to the background red-clump stars. This in- 
dicates that the cluster is located between these red-clump stars 
and the Sun. In other words, an upper limit to the distance can be 
established by comparing the target and control CMDs. As a re- 
sult, it can be concluded that the G6 1.48 +0.09 region is located 
at a distance less than ~4Kpc. It is worth mentioning that this 
rough estimation depends on the assumed intrinsic colors and 
magnitudes of a K2III star, as well as on the assumed extinction 
law. 

On the other hand, the local main sequence (MS) can be 
clearly seen in both CMDs, although for magnitudes fainter than 
Ks ~14, it appears more affected by extinction in the target field 
than in the control one. This additional extinction is due to the 
presence of the star forming region G6 1.48+0.09. At t his po int, a 
simplistic approximation can be done. iBinnev et al.l (l2000h find 
the age of the solar neighborhood to be 11 .2+0.75 Gyr If we as- 
sume that the solar neighborhood is composed mainly of stars 
with this age and solar metallicity, then a suitable isochrone 
can be over-plotted in Figure |6l In this context, a portion of a 
11.2 Gyr and solar metallicity isochrone, locate d at different dis- 
tance s has been over-plotted. The lAC-star ( A paricio & Gallar3 
l2004l) syn thetic CMD prog r am ha s been used to compute the 
isochrone. iPietrinferni et al] ( |2004|) stellar evolution library and 



ICastelli & Kuruc j ( 1200 3') bolometric correction library have 
been considered together with the Rieke et al. ( 1989) extinction 
law. Comparing these isochrones with the CMDs, it can be ob- 
served that the main difference between the target and control 
MSs starts at a distance of ~2.5 Kpc. That is, local MS stars lo- 
cated at a distance shorter than ~2.5 Kpc are not affected by the 
presence of the star forming region G61.48+0.09, but more dis- 
tant MS stars have an increased extinction when compared to the 
control field. 



As a result, this simple argument leads to a distance esti- 
mation of ~2.5 Kpc for the star forming region G61.48+0.09. 
In Section IH a more precise distance determination of 
G6 1.48 +0.09 is derived based on the spectroscopic analysis 
combined with photometric information, but it is interesting how 
this photometric analysis is able to provide a first approximation. 
Kinematical studies are compatible with distances of 2.0 and 6.5 
Kpc, so this photometric analysis points towards G61.48+0.09 
located at a distance closer to the smaller value. 

It is worth mentioning that no trace of a cluster MS is found 
in the G61.48+0.09 region's CMD. As it will be shown later, 
this cluster is affected by a high differential reddening, which 
explains the absence of a cluster MS. This differential redden- 
ing is shifting every cluster member in a different way along the 
reddening vector, so the cluster MS is highly blurred. 

3.2. Spectral classification 

The G6 1.48+0. 09 region's CMD has been used to select promis- 
ing cluster members, for which follow up spectroscopy has been 
done. The selection was done based on the location of these stars 
in the CMD and their position in the cluster In other words, the 
LS and MOS candidate list was performed based on photometric 
and spatial information. Note that this selection depends actually 
on the field morphology, as it was mainly observed with multi- 
object spectroscopy. 

The observed spectra are presented in Figure |7] Two ob- 
jects are classified as emission line objects, twelve turned 
out to be red objects, mainly G-type stars, and six resulted 
to be OB stars , from 09 to B3-B5. For the classi fi cation , 
the catalogs o f Hanso n & Contil (119961) . iHanson et all d 19981) . 
IWallace&Hhikle (.19971) ' and the relation between spectral type 
and eq uivalent width of the CO band head given in|p avies et aD 
("2007^ were used. The more recent atlas of Hanson et al.i (i2005l) 
or the high resolution atlas of Wallace & Hinkle ( 199g) were not 
used, as their resolutions (R~8000) are much higher than ours . 

Cool stars were first classified in spectral type using the 
[Wallace & Hinkld d 19971) catalog. To this aim, the CO band and 
the features around 2.265 fim (Ca/Sc/Ti), 2.207 (Na/Sc/V) and 
2.09 (CN) were used as primary indicators. Typical uncertainties 
(estimated from independent determinations by two of us) are 
+ 1 spectral subclass. After determining the spectral type, we de- 
termined t he luminosity class using the Spectral Type-EW(CO) 
relation in iDavies et al.l (|2007|) . This is illustrated in Figure |8] 
where we show the EW of the CO band head in our stars and 
Davies relationship. The CO EW was determined in the same 
way as indicated in iDavies et al. 1 (12007'): we measured the EW 
between 2.294 and 2.304 fim and adopted the average level be- 
tween 2.288 and 2.293 yum for the continuum. Errors were esti- 
mated from the local continuum rms. 

We see that the K and M stars fit clearly the giant calibration. 
The situation is less clear for the G stars, but their positions in the 
CMD indicate that they are not brighter than the blue stars in the 
sample. Therefore, they cannot be red-yellow supergiants. This 
same argument can be applied to stars #18 and #26, for which 
we could not get the spectra of the CO band (because of the 
usual displacement of the spectra when the slits are displaced 
from the mask central line), and for star #17 (F-type), where 
the CO-band is not visible. Additionally, we note that from the 
observed spectrum, star #17 could be a main sequence object, 
but then it would be too bright when corrected from reddening 
(except if it has a particularly strong differential obscuration). 
The spectra, magnitudes and CO band head equivalent widths 
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are thus compatible with all stars of type later than A, being 
yellow or red giants. 

Ea rly t ype spectra were clas sified using the lHanson & Contil 
(119961) and lHanson et al. l (fT998l) atlases. The main indicator used 
was the presence of Hel lines at 1.701 and 2.113. They are 
present in all stars we classified as ea rly-typ e st ars, except in 
star #2 2. According to^ Hanson & Contil (11996') and Hanson et alj 
(Il998l) these lines begin to disappear around B3/B5 in dwarfs. 
For supergiants the Hel 1 .70 line is present for all spectral types 
and the Hel 2.113 dissappears around B8/B9. The behavior of 
the Br series as compared with the other early-type stars points 
to a classification as B3/B5, probably of luminosity class V. Stars 
#4, #5 and #22 were thus classified as Bl, B2 and B3-B5 respec- 
tively. Star #12 shows in addition a trace of Hell 2.189 and was 
classified as 09. 

We tried to derive the luminosity class for the early types di- 
rectly from their spectra, but the low resolution in the H-band 
and the uncertainties in the reduction of the Br^ line in the K 
band prevented us from using the broadening of the H lines for 
the classification. We also tried to use the ratio of the Brl 1/HeI 
1.701//m. For star #12 this is difficult, as we are in the region of 
lower sensitivity for this indicator. Stars #4 and #5 have values 
favoring luminosity class V, while star #22 is clearly of luminos- 
ity class V according to this ratio. 

However, we can still use an indirect argument to assign a 
luminosity class to star #12. The stars #8 and #6 are red giants 
and they have very large extinctions. This indicates that they 
are behind the group of early-type stars. Therefore, they limit 
the possible distance to star #12 and, consequently, its absolute 
magnitude, which correspond then to a luminosity class V stars. 
The same argument can be applied to the rest of the stars in the 
cluster. Moreover, this conclusion is consistent with our previ- 
ous analysis of the CMD (see Sect. 13. It . If the early type stars 
were supergiants, then the derived cluster distance would be of 
the order of 7-8 kpc. This is not consistent with Fig.|6] where the 
CMD of the target field does not show the red clump strip for 
distances between 3 and 7 kpc. 

If the early-type stars were just giants the result would still 
be inconsistent. Star #22, and most probably also stars #4 and #5, 
are dwarfs. Assuming the stars #12 and #2 (that is classified as 
BOV below) to be giants would result in two different distances 
for the two groups. Stars #22, #4 and #5 are located at a distance 
of ~ 2.5 kpc, while stars #12 and #2 would be at 3.5 kpc. This 
argument leads us to a luminosity class V for all the early-type 
stars. 

The spectra of stars #1, #2, #3 and #14 show strong emission 
lines. The first three are in the central part of the cluster, and 
have been observed using a long slit that included all three stars. 
Unfortunately, the nebular emission in that region is very strong 
and highly spatially variable, so that at our low resolution their 
spectra are contaminated by emission lines, particularly in Br^. 
A similar problem arises with star #14, that was observed in the 
MOS mask and lies in the Bl subregion. We cannot separate the 
nebular and stellar lines. 

Star #3 displays a spectrum that resembles that of star #22 
(i.e., without any clear He features) except for the emission in 
the Bracket series. We classify it as a B3-B5 star contaminated 
by the strong nebular emission. 

Stars #2 shows a weak absorption at Hel 1.701 and a weak 
Brl 1 line, without trace of Hell absorption, which at this reso- 
lution is consistent with an 09 or BO spectral type. This is also 
consistent with the K-band spectrum showing no absorption fea- 
tures, with a mod erate emission i n Br^ and Hel 2.059, and with 
the analysis of iPuga et al.l (|2004|) that estimate a BOV spectral 



type for this star In our CMD this star falls between star #12 
and #4 that we classify as 09V and B IV. Therefore we classify 

star #2 as BOV. 

Star #1 (star #82 of iDeharveng et al.1 (|2000) and lPuga et J] 
(12004.) ) is the outstanding red object at the center of the cluster 
and it could be its main ionization source. In such case, its spec- 
trum would correspond to that of an early O dwarf. However, 
we find no trace of Hell lines or CIII/NIII/OIII emission at 
2.116, indicating that the spectral type is later than 07.5. At 
the same time, we see no CO absorption, so the star has to be 
of mid-F type or earlier If the extinction law towards the ob- 
ject is normal, then the spectrum and the observed K magnitude 
would be compatible with an early F supergiant, but then the 
star should be much redder in the color-color diagram. Therefore 
we conclude that the object has a spectral type between early-A 
and late-O. We cannot make a better classification, but we note 
that the K-band spectrum resem bles tha t of the YSO object in 
Gl 18.796-1-1.030 presented by Bik et al.l (|2005 ), so we classify 
it as an YSO. This is consistent with the strong reddening present 
in this star, and with the difficulties pointed out by Puga et al. to 
fit its SED. This classification, however, indicates that this is not 
the main ionization sourc e of the cluster, that is probably one 
of the hidden sources (see iPuga et aTll2004i for a more detailed 
discussion). 

Star #14 is quite similar to star #1 in the K-band, but with a 
stronger and broader emission in the Br^ line, as well as in the 
Bracket series in the H-band. We classify it as Be as we cannot 
identify further features. Stars #1 and #14 have very large (J- 
Ks) colors. For star #14 this is not a problem, because it still 
falls within our calibration stars range. Star #1 on the contrary 
falls beyond this limit, but we note that star #6, which falls in 
the same region, has infrared colors consistent with its spectral 
classification and a very large reddening. We thus think that the 
colors of star #1 are relia ble. The color excesses are cons istent 
with disks around star #14 dGehrz. Hackwell & Jones 119741) and 
star #1 ( Jiang et al. 20081) . 

The derived spectral types and some cross-correlations with 
previous works can be found in Table |2] 

4. Discussion 

The spectral classification of the early OB stars found, together 
with their photometry, allows the determination of various clus- 
ter parameters, such as its extinction, distance, age and mass. 



4. 1 . Extinction, distance and age 

Once the spectral type and luminosity class of a star is derived, 
its intrinsic color and absolute magnitude is known. Comparing 
them with its apparent color and m agnitude, the extin ction of this 
star can be derived. Assuming the Rieke et al. (1989) extinction 
law with R=3.09 (Rieke & Lebofskv.1985.) , the extinction, , 
of a star can be determined using the following equations: 



-H-Ks 



1.514 0.561 



(1) 



As the intrinsic colors of the cluster stars are known 
from their spectral classification, can be derived from the 
color-color diagram. Figure |9] shows the color-color diagrams 
for the stars in the G6 1.48 -1-0.09 target region. The solid and 
dotted lines show the reddening vector from the expected in- 
trinsic colors of a B 2V and G9III sta r, respectively, assuming 
the extinction law of iRieke et al.l (119891) . Filled and open circles 
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correspond to the stars classified in Section 13.21 as early (09V- 
B5V) and late (G3I11-M311I) type stars, respectively. It can be 
seen how the assumed extinction law reproduces very well the 
slope in the color-color diagram for both the early and late type 
stars. This suggest that the extinction law assumed for the dense 
molecular cloud is not very different from that in the interstellar 
medium. We also see that the early-type stars are affected by a 
strong differential reddening. The most reddened early-type star 
(#12) is actually far from what we would identify as the cen- 
ter of the cluster (centered at star #2). This is an indication that 
the cluster is still sweeping out the original molecular material. 
Besides, stars #l(YSO) and #14 (Be) present a high infrared ex- 
cess. 

We adopt the ICoxl (l2000h absolute V-magnitude and 
iDucati et al.l ( l200l i) IR colors for the early type stars except for 
star #12 (09V), for which Martins et al. (2005) observational 
scale colors have been adopted. The value of Afc^ has been ob- 
tained for each star using equation [T] and the mean of the two 
results is considered as the final A^^ . Table |3] lists A^^ results 
in column 1. It can be seen that the analyzed stars' extinctions 
vary from Ak^ = 0.9, or Ay ~ 8, to A/f, - 2.6, or Ay ~ 23. That 
is, the extinction is highly variable across the G61. 48-1-0.09 star 
forming region. 

Once the extinction of the individual stars is known, their 
distances can be derived comparing their apparent and absolute 
magnitudes. Column 2 in Table [3] lists the derived distances for 
the six cluster stars. The mean obtained distance is 2.5 Kpc. The 
uncertainty associated to the distance determination has several 
contributions. On one side, any uncertainty in the photometry of 
the stars, in their assumed intrinsic colors and absolute magni- 
tudes have an effect on the distance result. On the other side, the 
dominant source of error co mes from the extin ction. As it has 
been mentioned before, the feeke et il] (1 1989b extinction law 
with R=3.09 (Rieke & Lebofsky 198jrhas been adopted. The 
interstellar extinction l aw could be different in dense mo lecular 
clouds. For example, Breger, Gehrz, & Hackwell (Il98lb found 



an anomalous extinction law in Orion, and ICh inil (Il981b in the 
Ophiuchus dark cloud. So the interstellar extinction law for the 
cluster stars could be different from the adopted one, although 
not much according to Figure |9] If that is the case, any differ- 
ences between the adopted and real extinction laws would af- 
fect the derived distance. For this reason, we consider the sigma 
of these individual distances as the final uncertainty, that is, we 
adopt 2.5 + 0.4 Kpc as the distance of the cluster This distance, 
2.5 Kpc, is fully compatible with the shorter distance derived in 
previous kinematical studies and with the analysis of the CMD 
from previous section. 

Figure [TO] shows the calibrated CMD obtained for the star 
forming region G61.48+0.09. The spectroscopically observed 
stars are marked with filled circles (the numbers correspond to 
our own listing). In this CMD we can see the local main se- 
quence formed by slightly reddened field stars while cluster stars 
are expected to have larger reddening. This figure also shows the 
reddening coiTected position of the spectroscopically observed 
stars. Filled triangles show the position of the reddening cor- 
rected early type stars, that is, they show the reddening corrected 
cluster MS. Solid lines show the zero-age MS located at a dis- 
tance of 2.5 Kpc. Open triangles, on the other side, show the red- 
dening coiTected position of the late type stars. Vertical dashed 
lines show the color interval of the late type stars with spectral 
types between GOIII and M3III. 

Regarding cluster's age, only an upper limit can be estab- 
lished with present data. The earliest star found in this cluster is 
a 09V. If we assume that all stars in the cluster were formed at 



the same time, and we adopt t he calibration of stellar parameters 
given by Martins et al.'('2005') and the evolutionary models from 
Schaller et al. ( 1992), then we can conclude that the cluster must 
be younger than ~ 10 Myr. For older ages this 09V star would 
have evolved into a giant. This is consistent with the lack of red 
supergiants in the cluster, that nee d slightly less t han 10 Myr to 
develop from 25 M© stars in the ISchaller et al.l (1 1992ft models 
(including rotation in the evolutionary models does not change 
this conclusion). 

4.2. Mass 

We have estimated the stellar cluster mass of G61.48+0.09 as- 
suming a Salpeter IMF. Further, we also assume that the mass 
interval for 09V stars, the earliest type we hav e found, is 15.6 - 
18.8 Mq (from the observational scale in Ma rtins et al.ll2005h . 
With these assumptions the estimated mass contained in stars 
for G61. 48-1-0.09 is 1.5+0.5x10^ Mq, where the error has been 
estimated from our adopted uncertainty of +1 spectral subtype. 
This error is larger than the one derived from changing the spec- 
tral type-stellar mass relationship, but other important sources 
may contribute to enlarge the uncertainty and make our estima- 
tion a lower limit. 

The first one comes from the slope of the IMF. We have 
assumed a Salpeter form, which is in agreement with most 
determinations f or clusters in the mass range we investigate 
in G61.48+0.09 ( lMassevll2003h . Evidences of anomalou s flat- 
ter IMFs are found only in more massive clusters (iFigerl 
2005). Therefore we consider the Salpeter form adequate for 
G61.48+0.09. 

The second one comes from the possible presence of other 
early-type stars. We have seen in Section lTTI that the CMD stel- 
lar population above stars #2 and #12 is consistent with the 
nearby field population, which is also confirmed by the classi- 
fication of star #8 as a M2 III star. Therefore we don't expect 
a contribution of these stars to the early spectral type popula- 
tion of G61.48+0.09. More difficult is the possible contribution 
of star #1. If the strong IR excess is due to a disk around a 
young star, then its m ass is probably not higher than 20 Mq 
(ICesaroni et al.ll2006h . Assuming a second object of the same 
mass as star #12 we would obtain a cluster mass of 2900 M© 
(again with a Salpeter IMF). 

The final uncertainty would come from the lack o f an evident 
stellar ionizing source. The ionizing fluxes derived bv lPuga et al.l 
(|2004 cannot be fitted by stars #1 and #2 (Puga et al. stars 82 and 
83) at the spectral types and luminosity classes we have assigned 
them. Therefore, we have to look for another ionizing source, 
that presumably would be of earlier spectral type and higher 
mass. IPuga et al.l (l2004b propose that the source they identify in 
the L'-band as LI could be the source. Inspection of the Spitzer 
archive does not help, as the region is saturated at the long wave- 
lengths. Therefore, we emphasize that the mass we derive for the 
cluster is a lower limit. 



5. Conclusions 

A study of the stellar population of G61.48+0.09 has been done. 
This is part of our study of massive star cluster candidates in 
the Milky Way that is being carried out using LIRJS at WHT. 
New near-IR photometric data and, for the first time, near-IR 
spectroscopy of stars in the cluster region have been presented. 

From the study of the Color-Magnitude Diagram and the 
spectroscopic analysis, we conclude that G6 1.48 +0.09 is a 
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star forming region located at 2.5+0.4 Kpc, in agreement 
with authors that have chosen the short distance scale (see 
[peharveng et al. 2000). The cluster is younger than 10 Myr 
and we obtain a stellar mass of 1.5x10^ Mq, although with 
large uncertainty, as we cannot see the whole stellar content 
of G6 1.48+0. 09. In particular, our classification of observed 
spectral t ypes cannot account for the ionizing fluxes derived by 
iPuga et al . (2004) and it is unlikely that any of the stars observed 
in the near IR can be the ionizing source. Table|4]lists a summary 
of all cluster properties. 
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Fig. 1. Near- (Ks , LIRIS @WHT) and mid- ( 3.6 jum, GLIMPSE) infrared images of the star forming region G61 .48+0.09. The grey 
contours correspond to the lGarav et al.l (Il998l) continuum emission at 8.3 GHz. The positions of stars LI and L2 is also marked. 
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Fig. 2. Wide field DSS-2-infrared survey image of the G61.48+0.09 star forming region. The squares represent the observed target 
(a) and control (b) fields of view. 
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Fig. 3. 2MASS CMD of the target field. The stars used for the 
photometric caUbration are marked with open circles. 
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Fig. 4. Photometric calibration in the /-, H- and Ks -bands for 
the G61.48-I-0.09 field. Left panels show the photometric cali- 
bration from instrumental (INSTR) to the calibrated magnitudes 
(2MASS). Central panels show the difference between the cal- 
ibrated photometry and the 2MASS photometry. Right panels 
show the PSF fitting <t parameter provided by ALLFRAME. 
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Fig. 5. False color composition image of the central part of the 
G6 1.48 +0.09 field. The position of the spectroscopically ob- 
served stars si shown. The early-type stars (see Table |2]i have 
been marked with circles, and the rest with squares 
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(J-K3) 

Fig. 6. Color-Magnitude Diagrams of the target (a) and con- 
trol (b) fields. The black solid line shows the expected posi- 
tion of a K2III star as a function of its distance. To cons truct 
the K2III sequence, data for My h ave been taken from Coxl 
( I2OO O) and IR intrinsic colors from iDucati et al . (2001). Grey 
lines represent a fraction of a 11.2 Gyr and solar metallicity 
isochr one, located at different distan ces. The [P ietrinferni et al.l 
(|2004|) stellar evolution library and the lCastelli & Kurucz (2003^ 
bolometric correction li brary have been u sed to generate the 
isochron es. Besides, the [ Rieke et al.l (Il989l) extinction law and 
R=3.09 (lRieke&Lebofskvll985h have been assumed. 
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Fig. 7. H and K spectra of observed stars in the G61 .48+0.09 star forming region. 
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Fig. 8. EW of the CO band head versus spec tral type. Dot t ed line 
is the relation for supergiants presented by iDavies et al.l (l2007h 
and the dashed line is the corresponding relation for giants by 
the same authors. 
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Fig. 9. Color-color diagrams of the G6 1.48 +0.09 region (target 
field). Open and filled circles represents the late and early type 
stars, respectively, listed in Table |2l The solid line indicates the 
direction of the reddening v ector for a B2V star, while the dotted 
lines corre sponds to a G6III. Ducati et"an ( 1200 ll) IR intrinsic col- 
ors and the lRiekeetaII(ll989F extinction law have been adopted. 
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Fig. 10. Calibrated color-magnitude diagram of the G61.48+0.09 star forming region. The spectroscopically observed stars have 
been marked with filled circles. Filled triangles show the position of the reddening corrected early type stars, that is, they show the 
reddening corrected cluster main sequence. Open triangles show the reddening corrected position of the late type stars. The solid 
line shows the zero-age main sequence located at a distance of 2.5 Kpc. Vertical dashed lines show the intrinsic color interval of 
the late type stars with spectral types between GOIII and M3III. 
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Table 1. Imaging/spectroscopic observations details. 



Target 


RA (J2000) 


Dec (J2000) 


Date 


Filter/Grism 


Exp. Time (s) 


FWHM (") 


Standard star 


G61.48+0.09 


19 46 47 


+25 12 43 


July 22, 2006 


J Filter 


192.0 


0.6 


— 








July 22, 2006 


H Filter 


70.4 


0.6 










J Liiy zz, zuvju 


A^- FiiLCl 




n ^ 
u.z> 




Control field 


19 45 56 


+25 12 41 


June 6, 2007 


J Filter 


192.0 


0.9 










June 6, 2007 


H Filter 


70.4 


0.8 










June 6, 2007 


Ks Filter 


70.4 


0.9 




G61.48+0.09 


19 46 47 


+25 12 43 


June 5, 2007 


HK, R=945 (LS) 


1200 


1.0 


HD229700 








June 5, 2007 


K, R=2500 (LS) 


2400 


1.0 


HD229700 








Sept. 20, 2007 


HK, R=945 (MOS) 


2400 


0.7 


HD182761 








Sept. 20, 2007 


K, R=2500 (MOS) 


2400 


0.7 


HD182761 
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Table 2. Spectroscopically observed stars in G61.48+0.09. 



Star id 


RA (2000) 


(5(2000) 


J 


H 


Ks 


SpT 


Comments 


1 


19:46:47.585 


+25:12:45.64 


13.50 


12.06 


9.35 


YSO 


Ions-slit; #82 in Pusa et al. (2004) 


14 


19:46:47.823 


+25:12:29.98 


13.94 


12.34 


11.85 


Be 


12 


19:46:41.550 


+25:12:20.52 


14.60 


12.25 


10.92 


09 V 




2 


19:46:47.311 


+25:12:45.60 


12.60 


11.65 


9.47 


BO V 


Ions-slit: #83 in Puea et al. (2004) 


4 


19:46:48.993 


+25:13:28.86 


13.97 


12.46 


11.53 


Bl V 


5 


19:46:50.832 


+25:11:53.57 


13.06 


11.93 


11.25 


B2 V 




3 


19:46:47.033 


+25:12:45.55 


12.58 


12.52 


10.40 


B3-B5 V 


long-slit 


22 


19:46:44.779 


+25:12:46.70 


13.46 


12.54 


12.11 


B3-B5 V 




17 


19:46:50.608 


+25:13:11.86 


15.85 


13.93 


12.73 


F5 111 




9 


19:46:47.272 


+25:12:59.85 


15.96 


13.37 


11.77 


G3 111 




13 


19:46:47.118 


+25:12:34.05 


14.20 


12.80 


11.85 


G4111 




7 


19:46:49.037 


+25:13:00.72 


12.11 


11.44 


11.22 


G5 111 








1 T^. 1 1 . 1 O Q7 


lU.U / 


Q "iS 
y.jo 


Q 1 ^ 

y.Lj 


r^A TTT 
OO ill 




18 


19:46:51.209 


+25:14:11.36 


17.83 


13.76 


11.91 


G7111 




25 


19:46:40.728 


+25:10:57.91 


15.39 


13.24 


12.32 


G8111 




27 


19:46:46.924 


+25:12:13.37 


10.40 


9.64 


9.47 


G9111 




26 


19:46:40.380 


+25:13:21.48 


15.97 


13.10 


11.68 


KOlll 




15 


19:46:42.791 


+25:11:07.18 


14.06 


11.91 


10.96 


K3 111 




8 


19:46:44.697 


+25:12:13.08 


14.45 


11.00 


9.35 


M2 111 




6 


19:46:45.611 


+25:13:05.80 


13.77 


12.24 


10.74 


M3 111 
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Table 3. Derived extinctions and distances for the cluster stars. 



Star id 


Aks (mag) 


d (Kpc) 


12 


2.6 


2.1 


2 


1.7 


2.8 


4 


1.6 


3.1 


5 


1.3 


2.0 


3 


1.5 


2.1 


22 


0.9 


2.6 
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Table 4. Summary of cluster properties. 

Extinction Varies from Ak^ = 0.9 to 2.6 mag 
Distance 2.5±0.4 Kpc 

Age Younger than 10 Myr 

Minimum Mass 1.5x10' Mq 



